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We propose a method to study and constrain modified gravity theories for dark matter using 
CMB temperature anisotropics and polarization. We assume that the theories considered here have 
already passed the matter power-spectrum test of large-scale structure. With this requirement met, 
we show that a modified gravity theory can be specified by parametrizing the time evolution of its 
dark-matter density contrast, which is completely controlled by the dark matter stress history. We 
calculate how the stress history with a given parametrization affects the CMB observables, and a 
qualitative discussion of the physical effects involved is supplemented with numerical examples. It 
is found that, in general, alternative gravity theories can be efficiently constrained by the CMB 
temperature and polarization spectra. There exist, however, special cases where modified gravity 
cannot be distinguished from the CDM model even by using both CMB and matter power spectrum 
observations, nor can they be efficiently restricted by other observables in perturbed cosmologies. 
Our results show how the stress properties of dark matter, which determine the evolutions of both 
density perturbations and the gravitational potential, can be effectively investigated using just the 
general conservation equations and without assuming any specific theoretical gravitational theory 
within a wide class. 

PACS numbers: 04.50.Kd, 95.35.-|-d, 98.70.Vc 



I. INTRODUCTION 

Various cosmological observations indicate that almost 
95% of the energy content of our Universe is in some dark 
form that can only be detected through its gravitational 
effects. The standard explanation for this "dark sector" 
involves some kind of cold dark matter (CDM), such as 
heavy weakly interacting particles or primordial black 
holes, and dark energy, which could be a manifestation of 
a cosmological constant (A) or some exotic matter field. 
This standard hybrid picture (ACDM) has so far passed 
several different cosmological tests, and is dubbed the 
"concordance model". Despite of its successes, however, 
it is not without problems in accounting for structure on 
galactic scales - for example, the greatest challenge is 
to form bulge-free bright spiral galaxies and dwarf spiral 
galaxies. The galaxies in CDM have either a dominating 
stellar bulge or a dominating CDM cusp at the center. 
In essence LCDM over-predicts the dark-matter effects 
required by the empirical TuUy-Fisher relation of spiral 
galaxies - and is in need of an experimental identification 
of the CDM particles with an understanding why the 
dark energy possesses its particular energy density. 

Since dark energy affects cosmological models mainly 
via gravitational effects, it is possible to imagine that the 
effects of an explicit material source for the dark energy 
can be mimicked by a change in the behaviour of the 
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gravitational field at late times. This is an unusual re- 
quirement since we have expected deviations from general 
relativity to arise in the high spacetime curvature limit 
at early times rather than in the low spacetime curvature 
limit at late times. However, we should note that the ad- 
dition of an explicit cosmological constant to general rela- 
tivity, as in the concordance model of ACDM, is already a 
particular example of such a low spacetime curvature cor- 
rection. There have been many investigations of gravita- 
tional alternatives to dark energy pL, 2, 3, 4, 5, 6, 7, 8, !|, 
and it appears that most of these attempts create dif- 
ferent problems in either local gravitational systems 
or large-scale structure (LSS) formation, or even both 
ll PJIl, [11, [13, [III, [11 [11 111, 111 
20, [Sll on possible ways to over- 
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come these problems). In response to these investiga- 
tions, frameworks to test modified gravitational dark en- 
ergy models have also been developed [32, 33, 34, 35, 3(|. 

The situation is nonetheless different in the dark mat- 
ter arena, where the leading modified gravity model, Mil- 
grom's modified Newtonian dynamics (MOND) [33], ap- 
peared more than two decades ago, but lacked a con- 
vincing general-relativistic formulation with a set of rel- 
ativistic field equations that are applicable to cosmol- 
ogy. This hurdle has recently been overcome by Beken- 
stein's tensor-vector-scalar (TeVeS) model which re- 
duces to MOND in the relevant limit. Actually, what 
the TeVeS model provides is more than a relativistic 
framework to investigate cosmology - it has been found 
that in TeVeS the formation of LSS, which was thought 
of as a problem for modified gravity theories before, 
could also be made consistent with observations ^3^, i4^ . 
Meanwhile, TeVeS has also been shown to work well on 
smaller scales {e.g., mimicking cold dark matter in strong 
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gravitational lensing systems [4l|, \^ , producing ellipti- 
cal galaxies, barred spiral galaxies and even tidal dwarf 
galaxies [13, ISl , being consistent with solar system tests 
and so on) and inherit the advantages of MOND over 
the CDM model 45] {e.g., in explaining the TuUy-Fisher 
law and galaxy rotation curves), this discovery attracts 
much interest on TeVeS and its generalizations. Also, 
since TeVeS manages to grow large-scale structure us- 
ing the growing mode of its vector field, this stimulates 
the investigation of vector-field cosmology [4611 in general 
dllis, 49, 50,[5ll,[52,[5l,[5i,[5ij (in Ref. for exam- 
ple, a very general vector-scalar field framework, the so- 
called dark fluid, is proposed which can reduce to many 
existing models in appropriate limits), and it is found 
that the LSS in these theories could also be consistent 
with observations [56]. 

Despite of this encouraging success, one should bear 
in mind that LSS only provides one test of any structure 
formation model. In fact, the matter (or galaxy) power 
spectrum we observe today only reflects the large-scale 
matter-density perturbation today (Smo), rather than 
its evolutionary history. Thus, even though the matter 
power spectrum, P{k), is compatible with observations, 
the evolution path of Sm may well be different from that 
followed in the ACDM paradigm. This situation is shown 
particularly clearly in Fig. 3 (lower panel) of po| . where 
the growth rate of baryon density perturbations is en- 
hanced only at late times. The different evolution history 
of matter-density perturbations may have significant im- 
pacts on various other cosmological observables, such as 
the cosmic microwave background (CMB) power spec- 
trum, and influence the growth of nonlinear structure. 
We would like to determine what these impacts are. 

In this paper we take a first step in that direction. 
We will concentrate on the influences of general mod- 
ified gravity dark-matter models on CMB observables. 
Our main assumption is that the modification to gen- 
eral relativity (GR) can be expressed as an effective dark 
matter (EDM) term and moved to the right-hand side of 
the field equations so that the left-hand side is the same 
as in GR. This EDM term, hke the standard CDM, gov- 
erns both the background and perturbation evolutions. 
This assumption is justified by observing the fact that 
in TeVeS, as well as in the general vector-field (or f{K)) 
theories, the terms involving the vector field are essen- 
tially just the EDM term described here. Furthermore, 
we make some simplifying assumptions. First, any ex- 
plicit dark energy is neglected. The main effects of dark 
energy (if not too exotic in origin) are to modify the back- 
ground evolution at late times and cause the decay of the 
large-scale gravitational potential. Neglecting it does not 
affect the essential features of our model but will sim- 
plify the numerics greatly. Second, the background evo- 
lution is exactly the same as that in the standard CDM 
(SCDM) model, which should also be a good approxima- 
tion. Third, the matter power spectrum, or equivalently 
5mO, is the same as that of SCDM, because any deviation 
from the latter should be stringently constrained by LSS 



observations, and because both TeVeS and f{K) mod- 
els have claimed to reproduce the observed matter power 
spectrum. Hence, we are fixing 5„m and investigating 
how different evolutions of Smio.) (where a is the cosmic 
scale factor) affect the CMB power spectrum. Our the- 
oretical framework is designed to be more general than 
is required for this purpose alone, and could be used to 
investigate features of other cosmological models. 

This paper is organized as follows: in § |TT] we set out 
the theoretical framework for our investigation and in- 
troduce more details of the cosmological model. In § IIIII 
we describe briefly how a general dark-matter component 
affects the CMB power spectrum and then, in § IIVI sup- 
plement this discussion with a numerical example. We 
consider three special cases which span all the possibili- 
ties in the model, and explain them one by one. Finally, 
in §|V]we provide a summary of our results together with 
some further discussion of the assumptions on which they 
are based. 



II. THE THEORETICAL FRAMEWORK 

In this analysis we use the perturbed Einstein equa- 
tions in the covariant and gauge invariant (CGI) formal- 
ism. 



A. The Perturbation Equations in CGI Formalism 

The CGI perturbation equations in general iE-theories 
are derived in this section using the method of 3 + 1 
decomposition ^5?|. First, we briefly review the main 
ingredients of 3-1-1 decomposition and their application to 
standard general relativity [s^l for ease of later reference. 

The main idea of 3-1- 1 decomposition is to make space- 
time splits of physical quantities with respect to the 4- 
velocity u° of an observer. The projection tensor hab is 
defined as hat = Qab — UaUb and can be used to obtain 
covariant tensors perpendicular to u. For example, the 
covariant spatial derivative V of a tensor field Tj^'..^ is 
defined as 

v'^Tl.i^h'ih)--- hih^,--- hy^T^:::^. (i) 

The energy-momentum tensor and covariant derivative 
of the 4- velocity are decomposed respectively as 

Tab = T^ab + 2q(aUb) + pUaUb - phab, (2) 
Valifc = (Tab + -^ab + ^Ohab + UaAb. (3) 

In the above, nab is the projected symmetric trace-free 
(PSTF) anisotropic stress, qa the vector heat fiux vec- 
tor, p the isotropic pressure, aab the PSTF shear ten- 
sor, zuab = V^qU;,], the vorticity, 9 = V^Uc = ia/a (a is 
the mean expansion scale factor) the expansion scalar, 
and Ab = Ub the acceleration; the overdot denotes time 
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derivative expressed as 4> = u°'V a(j}, brackets mean an- 
tisymmetrisation, and parentheses symmetrization. The 
4- velocity normaUzation is chosen to be u°-Ua = 1. The 
quantities 'Kab,Qa, P,P are referred to as dynamical quan- 
tities and Uabi'^^abTOT^a as kinematical quantities. Note 
that the dynamical quantities can be obtained from the 
energy-momentum tensor Tab through the relations 

p - TabU''u\ 
P = -W'Tab, 



T^ab = Kh'^Tcd + phab- 



(4) 



Decomposing the Riemann tensor and making use the 
Einstein equations, we obtain, after linearization, five 
constraint equations [s^ : 
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W^b + -6qa 
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and five propagation equations, 
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(7) 
(8) 
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Here, Sabcd is the covariant permutation tensor, £ab and 
Bab are respectively the electric and magnetic parts of 
the Weyl tensor Wabcdi defined by £ab = u'^u'^Wacbd and 
Bab = —^u'^u'^'Sac^yVefbd- The angle bracket means tak- 
ing the trace-free part of a quantity. 

Besides the above equations, it is useful to express the 
projected Ricci scalar R into the hypersurfaces orthogo- 
nal to M° as 



R 



^ 3 



(15) 



The spatial derivative of the projected Ricci scalar, yy^ 
^aVa-R, is then given as 



2a - 

Va = nVaP - —9VaO, 



(16) 



and its propagation equation by 
2 
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Va + y?7a 



6'aVaV • A - okVqV • q. (17) 



Finally, there are the conservation equations for the 
energy-momentum tensor: 

p + (p + p)e + v\ia = 0, (18) 

qa + ^eqa + {p+p)Aa-VaP + y''T^ab = 0. (19) 

As we are considering a spatially-flat universe, the spa- 
tial curvature must vanish on large scales and so i? = 0. 
Thus, from Eq. (fT5|l . we obtain 
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np. 



(20) 



This is the Friedmann equation in standard general rel- 
ativity, and the other background equations (the Ray- 
chaudhuri equation and the energy-conservation equa- 
tion) can be obtained by taking the zero-order parts of 
Eqs. (Unmil), yielding: 



'+^^' + f(p + 3p) - 0, 
p+{p + p)0 ^ 0. 



(21) 
(22) 



All through this paper wc only consider scalar pertur- 
bation modes, for which the vorticity Wab and magnetic 
part of Weyl tensor Bab are at most of second order [s^l , 
and so are neglected in our first-order analysis. 



B. Perturbation Equations in fc-space 

For the perturbation analysis it is more convenient to 
work in the k space because we confine ourselves in the 
linear regime and different fc-modes decouple. Following 
[s?! , we shall make the following harmonic expansions of 
our perturbation variables 

VaP=Y.-XQl VaP=J2-^'Qa 

^ — ' a — ' a 

k k 

k k 

VJ = J2'^ZQ'a '^ab-Y.-'^Q'b 

k k 

i/aa = V khQl Aa=y^ -AQ>1 

k k 
£ab^-Y.K^Qab Va = J2'i^Qa (23) 



in which Q*^ is the eigenfunction of the comoving spatial 
Laplacian a^V^ satisfying 
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and g^,Q^, are given by = fVaO^Q^, = fV^.Q^^ 
In terms of the above harmonic expansion coefficients, 
Eqs. (0 11 [m [131 [m [ni) can be rewritten as [13] 



Kqa 



-KO^ [k{\l + X)+mq] 



k{a' +n(7) = k^{(t) + A) - -na^n, 



(24) 
(25) 
(26) 



k\<j)' + H(I)) 



-Ka^ [k{p + p)a + kq - W - HU] , (27) 



k^rj = KXa^ - 2kHZ, 
kri' = -Kqa^ - 2knA 



(28) 
(29) 



where TC = a' /a = ^aO and a prime denotes the deriva- 
tive with respect to the conformal time r [adr = dt). 
Also, Eq. (fT9|l and the spatial derivative of Eq. p8)) be- 
come 



q' + AUq +{p + p)kA - kX^ + -fcH 

3 

X' + 3h'{p + p)+ m{X + XP) + kq 



C. The Main Equations 



0, 
0. 



(30) 
(31) 



Recall that we are treating the modifications to GR as 
an EDM term and so include them in the (generalized) 
energy-momentum tensor Tat to maintain the standard 
form of the Einstein equations. Thus, we can distinguish 
its different components pedm,Pedm, ga,EDM, 7ra6,EDM 
and their conservation. Depending on the specific model, 
the expressions for these components can be very differ- 
ent, but their conservation equations will take the same 
form. In particular, since the EDM has no coupling with 
standard model particles such as photons and baryons, 
it satisfies a separated conservation equation, Eqs. pO[ 
ED: 



v'T^UM+'KvEBM + kA-kAp + -kA^ = 0, (32) 



A' 



EDM 



kZ - 3HA + 'SHA„ + kir 



EDM 



= (33) 



where we have defined the EDM peculiar velocity wedm = 
Qedm/pedm, the density contrast Aedm = -^edm/pedm, 
and Ap = X^-^-^/peum, A^ = Hedm/pedm (cf. ^ IIIB|) 
for the EDM, and used the fact that pedm = to re- 
produce the standard CDM background evolution. The 



prime here is the derivative with respect to the conformal 
time and Ti. — a' /a. 

On the other hand, the spatial derivative of Eq. pO|) 
gives the evolution equation for Z as 



k'Z + knZ - kA+^K{X + 3X^)0^ =0 (34) 



and in this equation X^ X^ are respectively the density 
and pressure perturbations for all the matter species, in- 
cluding the EDM. For convenience we shall work in the 
frame where yl = 0. In this case, if the universe is dom- 
inated by the EDM, then the three equations above can 
be combined to eliminate Z and uedm^ 



A" 

^edm 



HA 



1 2 

EDM - -^KpEDMO- Aedm - -^k A^ 



-37^A' 



A„ 



0(35) 



If other matter species could not be neglected, as is in the 
radiation dominated era and early matter era, then we 
only need to correct the above equation by adding to it 
^'tpEDMa^AsDM another term with ^Kpi,a^ Ai,+ nprO^ Ar 
(where pf,, p^. are the baryon and radiation energy densi- 
ties) which comes from the last term in the left-hand side 
of Eq. ([Ml and these new terms do not affect the qual- 
itative features of our discussion (we will include them 
in the mmierical calculation). The Eq. ([55[) tells us that 
the evolution of the density perturbation in EDM is com- 
pletely controlled by its stress history: the EDM equa- 
tion of state (EOS) wedm = Pedm/pedm controls the 
background expansion and is set to zero here; the other 
two stress variables Ap, Ajr then drive the evolution of 
Aedm through Eq. ([35[) . Note that the stresses of the 
EDM are external functions determined by microphysics 
(for particle dark matter) or a particular modified grav- 
ity theory, and must be specified to close the system of 
Einstein equations and conservation equations. In the 
special case where Ap = A^r = 0, we reduce to the CDM 
model for which Acdm oc a in the matter era; but in gen- 
eral there will be deviations from this growing solution. 

Next, we look at the evolution of the gravitational po- 
tential, (see ? [ira[ITB[) . By manipulating Eqs. P51 [^ 
Wf\ [5D[) . and working again in the ^ = frame, we can 
eliminate the terms involving q and obtain the following 
evolution equation 
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2n' + [1 + 3 



-^pa ( Ap - -A ) - — «;pa A, - ^t^pa 



2fc2 



Kpa 



5 + 3^ 



7^ + - I H 

P. 



P' 



a: 



A^ 



(36) 



where p, p are the total energy density and total pressure 
respectively and p = pedm, with p = for an EDM dom- 
inated Universe. The term ^upa^ ^Ap — ^A^ is zero for 

radiation and baryons, but in general could be nonzero 
for the EDM. Again, we see that the stress history of 
the EDM completely specifies the evolution of cj). As we 
shall see in what follows, the dark matter component 
influences the CMB power spectrum through the gravi- 
tational potential, and so the stress history of the EDM 
is of crucial importance for the CMB. 

In summary, we have seen that the properties of the 
(isotropic and anisotropic) stresses of the EDM determine 
the evolution of the matter density perturbation and the 
gravitational potential, and thereby determine the pre- 
dicted matter power spectrum (through the former) and 
the CMB spectrum (through the latter) |58|. As we dis- 
cussed in § m the predicted matter power spectrum of a 
theory records the matter energy density perturbation at 
a specific time - today, and it is easier to make it con- 
sistent with observations, as evident from the studies of 
TeVeS and f{K) theories. Thus, in this work we shall 
start from an EDM theory which is constructed already 
to predict an acceptable matter power spectrum; that is, 
we fix the matter energy density perturbation today and 
calculate the influence of different Sm evolution paths on 
the CMB spectrum. In this way we can reduce our model 
space to one that is of realistic interests. 



III. THE DARK MATTER EFFECTS ON CMB 

It is helpful to have a brief review of the CMB physics 
and how it is affected by the dark matter before we go into 
the numerics to show the effects of EDM stresses on the 
CMB. Here we will just present a minimal description of 
this topic, for more details see the reviews [59l[60ll6ll[6^ . 

The primary CMB spectrum is determined by inhomo- 
geneities in the CMB photon temperature at the time of 
recombination. Prior to the recombination, photons cou- 
ple tightly by Thomson scattering with electrons which 
themselves couple to baryons via Coulomb interaction, 
thus to a first approximation photons and baryons com- 
bine into a single fluid. Dark matter, on the other hand, 
does not couple electromagnetically, but only interacts 
through gravitational effects and so contributes to the 
gravitational potential in which the baryon-photon fluid 



moves. 

When the density perturbation in the photons grows, 
hot (cold) spots appear where the local photon density 
is higher (lower) than average. The higher photon den- 
sity means higher pressure, which will tend to counteract 
the growth of local photon density. This will lead to 
acoustic oscillations of the gauge invariant photon tem- 
perature perturbation O, which is also a characterization 
of the photon number density perturbation. In reality 
this picture becomes a bit more complicated because of 
the interplay between baryons and dark matter: baryons 
couple tightly to photons, and move with them so that 
the inertia of the fluid is increased, while the gravitational 
potential produced by dark matter drives the oscillations. 
More specifically, if we use the multipole decomposition 

e 

where Pi is the Legendre function and fc/i = k • 7 with 
7 being the direction of photon momentum, then in the 
tight-coupling limit the monopole Qq satisfies a driven- 
oscillator equation j63j : 

©0 + Y^WQo + fc'c^eo = F, (37) 

in which 

F = ^n<S>'-lkH, (38) 

where 

i? EE 3ph/4p^, = 1/3(1 + i?) 

and 

1 1 a2 

are, respectively, the (frame-independent) Newtonian po- 
tential and curvature perturbations. The dipole moment 
Oi, which equals the peculiar velocity of baryons thanks 
to the tight coupling, satisfies fc8i -I- 3HQq + 3$' — 0, 
and all higher moments Qe {£ > 2) vanish because the fre- 
quent scattering makes the photon distribution isotropic 
in the electron rest frame. 

Thus, on large scales the photon temperature pertur- 
bation displays a pattern of driven and damped oscilla- 
tion. On scales smaller than the photon mean free path. 
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which itself grows in time, however, the tight coupHng 
approximation is no longer perfect and quadrupole mo- 
ments of the temperature perturbation needs to be taken 
into account. This introduces a dissipation term into the 
oscillator equation above, which damps the oscillations. 

At the time of recombination, the number density of 
free electrons drops suddenly and there ceases to be cou- 
pling between baryons and photons. The CMB photons 
then free stream towards us. This free-streaming solution 
is given by [1^ 

e,(ro) « (eo + «')(n)(2^-^l)j,(fcAr,) 

+ei(T*) [^j,_i(fcAn) - (£+ l)j,+i(fcAT,)] 

+{2e + 1) J^\^' - ^')n[k{To - T)]dT (39) 

where tq, are the conformal times today and at recom- 
bination, Ar* = To — T* and je{x) is the spherical Bessel 
function of order £. The terms -t- 4" and Qi are re- 
spectively the monopole and dipole moments of the CMB 
temperature field at , and can be calculated using the 
equations above (a ^E" is added to 6o, which accounts for 
the redshift in the photon energy, and thus temperature, 
when it climbs out of the Newtonian potential). Finally, 
the CMB temperature-temperature spectrum is defined 
as 

J 21+1 (40) 

These are the leading-order effects in the CMB physics, 
and we could see how they affect the CMB power spec- 
trum. Because jeix) peaks strongly at ^ ~ a;, so pertur- 
bations on very large scales (small k) mainly affect the 
low-^ spectrum. If the Universe is completely dominated 
by matter at then the first term in Eq. p9p is given 
by (Oq -I- ^')(t*) « ^'(r*)/3, and accounts for the ordi- 
nary Sachs- Wolfe (SW) effect; meanwhile, if the potential 
^ — ^ = 2(j> decays between and tq, then the integra- 
tion in Eq. p9p will also make a significant contribution 
as photons travel in and out of many time-dependent po- 
tentials along the line of sight, and this is the integrated 
Sachs- Wolfe (ISW) effect. 

Going to higher i one can see the peak structures of 
the CMB spectrum. The peaks appear since at r* , when 
the oscillating pattern of 6o [c.f. Eq. ([57)1 ] freezes, 6o 
might be just at its extrema for some scales {k); and these 
extrema are then converted to extrema of Ce through 
Eqs. SOI)- Since what appears in Eq. gO]) is \Qe\'^, 
both the maxima and minima of (0o + ^')(t») will appear 
as peaks in Cg. If there are no baryons and the potential 
^I^ is constant, then the even and odd peaks should be of 
the same amplitude. The inclusion of baryons effectively 
increases the inertia of the baryon-photon fluid and dis- 
places the balance point of Gq -I- ^, and as a result, after 
taking | • • • P, the odd and even peaks appear to have dif- 
ferent heights. At the same time, from Eqs. ([551 HO)) , we 
see that the dipole 0i(t,) also contributes to Ci through 
modulation. But Eq. p9p shows that its power is more 



broadly distributed and so its contribution is significantly 
smaller than that of the monopole. Where the monopole 
vanishes the dipole becomes important and this is why 
the troughs of Ci are of nonzero amplitude. 

If one goes to still higher £, there are two effects. First, 
the higher i moments mainly receive contributions from 
the small-scale (large k) perturbations, which began os- 
cillating already during the radiation-dominated era (not 
much earlier than ) . Since in the radiation era the po- 
tential <f> decays, this decay will drive the monopole oscil- 
lation through the source term in Eq. ([57)1 and lead to an 
increase in the oscillation amplitudes. This explains why 
in some models the third peak is higher than the second. 
Second, as discussed above, on very small scales there 
is severe damping of the oscillations because of photon 
diffusion. The combination of these two effects causes 
the CMB power in is higher than the third peak to be 
significantly damped. 

We can now highlight the key places where the dark- 
matter effects enter CMB physics through the gravita- 
tional potential it produces and contrast the situation 
with that when EDM is employed as a substitute. First, 
if CDM is replaced by EDM, then in the analysis of § [II] 
the time evolutions of 0, and thus of and ^' are mod- 
ified. If the modification only becomes significant after 
r,, the ISW effect could be different from standard CDM 
(where it effectively vanishes) . If it differs from standard 
CDM before then the SW effect will be altered as 
well. Second, the deviation of <& and 5" from their values 
in standard CDM before t* can change the zero point of 
the monopole oscillation through baryon loading, thereby 
modifying the relative heights of the odd and even peaks. 
Third, the modification of $ and 5* also changes the driv- 
ing force in Eq. (j37p . and results in different amplitudes 
for the at high £ . In § IIVI we shall give some numeri- 
cal examples showing these effects explicitly in our EDM 
model. 



IV. NUMERICAL EXAMPLES 

In this section we turn to some numerical examples to 
illustrate the qualitative analysis of § IIIII As mentioned 
in § n we choose to fix the endpoint of the evolution 
of the matter energy density perturbation. This gives 
us a freedom to parametrize the evolution of AEDM(a), 
which is generally different from that of Acdm (a) oc a (in 
matter regime) . Once the Aedm parametrization is given 
(or in other words the model is set up) , Eq. (|35p becomes 
an evolution equation for Ap and Att . This is insufficient 
to solve for Ap and A^, so in the following numerical 
calculations we consider three cases: (1) A^ — 0, (2) 
Ap = 0, and (3) a more reahstic case where both Ap and 
Att are nonzero. Since any scale dependence of Aedm will 
alter the shape of the matter power spectrum, and thus 
incur stringent constraints, we parametrize Aedm (a) so 
that it is independent of k. 
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A. The Case of = 



When Ajr = and Aedm is specified, Eq. (|^ becomes 
a first-order evolution equation for Ap 



3HA' 



37i' + 37Y^ - -KPEDMfl^ 



5, (41) 



where the source term S is given by ^ 



S — — AgQM — "^^EDM + 2^^EDMa^AEDM 
+ — Kpba^ At + KprO^Ar- 



(42) 



The quantity Ap does not appear directly in the expres- 
sion for gravitational potential [c.f. Eqs. ([25l [27|) ]. how- 
ever it does affect the potential (f) indirectly through the 
evolution of vebm Eq. ((32)) . Consequently, there are now 
two new variables (Ap and ■^edm) to evolve in the model 
and we need to specify their initial conditions. 

Here, we adopt the simplest and most direct approach, 
namely to assume that the EDM evolves as CDM prior 
to some initial time r^, when the scale factor is a.; and 
the dark matter density perturbation A^, and for a > a; 
the dark matter density perturbation begins to evolve as 
Aedm (a)- This means that for a < the variables Ap 
and wedm ^ will remain zero: this provides the desired 
initial conditions. This simple approach captures most 
of the interesting features about EDM. However, by as- 
suming that at early times EDM is just CDM, we cannot 
account for more complicated issues such as the primor- 
dial power spectrum, and we will comment on this in the 
concluding section. 

Now we could turn to a specific example of Aedm (a) 
parametrization. We let Aedm (a) equal AcDM(a) in the 
standard CDM model at times a; and af, and assume 
that the deviation from standard CDM occurs between 
and a/. Here, a/ < ao and ag = 1 is the current time. In 
order to characterize the deviation from standard CDM 
between at and af, we introduce a parameter b to denote 
the ratio of Aedm to Aqdm at the mean time a = (ui + 
a/)/2. Since for standard CDM we have AcDM(a) « 
a, the parametrized Aedm is simply a parabola which 
passes through three points (a.;, A^), (a/, AiUf/ai) and 
/Oj+o/ ^ Ai a,+aj |-,g^.^yggj^ times when the scale factor 
lies between Ui and a/: 



AEDM(a) = — Q-2— , (b-l)(a-a,)(a-a/) 

(43) 

for a G [fli, af]. Some examples of Aedm (a) with differ- 
ent choices of h are shown in Fig.[TJ It clearly shows that 




0.2 0.4 0.6 

Scale factor a 



FIG. 1; Some example parameterizations for the evolution of 
Aedm, normalized to its current value Aedmo, as described 
in Eq. (|43p . The values of the parameter 6 are indicated; the 
other parameters are = 0.005 and aj = 1. 



b characterizes the deviation from a standard CDM evo- 
lution (5 = 0). Note that the parametrization Eq. ([43|l 
is just an example for illustration, and specific modified 
gravity models like TeVeS and f{K) theory may lead to 
different parametrization which however can similarly be 
tested. 

Then with the parametrization of AEDM(fl) in Eq. (j43p . 
we could numerically evolve the relevant perturbation 
equations and see the changes in the CMB power spec- 
trum. First, we assume > t*. In this situation the 
CMB physics prior to last scattering is unaffected and so 
from Eq. ([39]) we see that only the integration is mod- 
ified. This is because in the standard CDM model the 
gravitational potential, </) = (^f — ^)/2, remains constant 
during the matter epoch, and the integrand in Eq. 
vanishes. For the present EDM case, however, Eq. 
dictates that 4>' ^ even in the matter era; consequently, 
there will be a significant ISW effect to boost the low-£ 
CMB power. As shown in Fig. [21 the earlier the EDM 
evolution deviates from that of CDM {i.e., the smaller 
ai is), the earlier (j) begins to evolve and the larger the 
cumulative ISW effect is. 

Next, consider the case where the deviation from CDM 
evolution begins earlier than last scattering. In this 



^ Here we include the effects of baryons and radiation explicitly. 
^ Recall that we work in the A = frame, in which t^cDM = 0- 



^ Note that in reality the universe is not completely matter dom- 
inated between a = 0.0002 and a = 0.002, and as a result the 
parametrization Eq. 1 143 II . with fe = 1, is not exactly the same as 
Acdm(i). The qualitative features, however, are not affected by 
this small deviation. 
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FIG. 2: (color online) The CMB power spectrum for our Att ~ 
model with AEDM(a) parameterized as in Eq. (|43|l . Three 
values of b are adopted: b — 1.2 (blue curves), 6 — 1.1 (red 
curves) and & = 1.0 (black curve). The dashed curves are the 
cases where ai = 0.005 and the solid coloured ones = 0.05. 
For all the curves we set a/ = 1.0. 



FIG. 3: (color online) The primary CMB spectrum with the 
ISW contribution removed, for the parameters ai — 0.0002, 
a/ — 0.002 and 6 — 1.2 (blue curve), 1.0 (black curve) and 
0.8 (red curve) respectively. 



case the evolution of Aedm and thus the potential is 
changed before last scattering, and correspondingly the 
first two terms in Eq. ([55)1 are modified as well. In Fig. [3] 
we have displayed the primary CMB power spectrum for 
the Aedm (a) parametrization Eq. (j43p with parameters 
a, = 0.0002, a/ = 0.002 ^ ,6=1.2, 1.0, 0.8 and the ISW 
effect switched off for simplicity. A larger value of b im- 
plies a larger EDM density perturbation (for the scales 
relevant to the first acoustic peak) before the last scatter- 
ing and therefore a deeper Newtonian potential "if. This 
means that the CMB photons experience larger redshifts 
when climbing out of the potential after last scattering, 
and the effective temperature (6o + ^) ('''*) is lower, lead- 
ing to a suppressed first CMB acoustic peak. When b is 
smaller the opposite effect occurs (c.f. Fig. [3]). At the 
higher-^ peaks this effect is less significant because the 
potential for the relevant (smaller) scales has already de- 
cayed. However, as discussed in § IIIII the change the 
potential prior to last scattering also modifies the equi- 
librium point of the monopole oscillation and alters the 
relative heights of odd and even CMB peaks. If ^ is 
a constant, the difference between odd and even peaks 
in |8o + 5*1 is proportional to \'^\ [60], so increasing b 
amplifies this difference by increasing 5", and as shown 
in Fig. O the third peak becomes higher and the second 
peak becomes lower for larger values of b and go oppo- 
sitely for smaller b. In this figure the fourth peak and 
onwards are not affected significantly because our param- 
eters are chosen conservatively. There however could be 
alternative dark matter models where the higher peaks 



also deviate from the CDM predictions, for an example 
see Fig. 4 (upper panel) of |39| . 

The CMB power spectrum has been measured to high 
precision by many experiments (see for example (b^. [65j) 
and will be further improved in the future, so it could 
be used to constrain the EDM model here. For higher £s 
the CMB data could be used directly. For lower £s its us- 
ability is limited by the cosmic variance. If the deviation 
is significant (such as those in Fig. [2]) , then we could use 
CMB data alone to constrain the model, as in [6^ . Oth- 
erwise, we could cross correlate the observed ISW with 
the matter density perturbation observable [ij ; because 
both are modified in the EDM model, we should expect 
different correlations from those created by CDM. Actu- 
ally, this technique has been applied to TeVeS [1^ but 
may not serve as a perfect EDM discriminator because 
there will be contaminations from the dark energy com- 
ponent in the galaxy-CMB correlation. Meanwhile, the 
general EDM model has different gravitational potential, 
matter density perturbation, as well as a possible dif- 
ferent redshift distribution of lensing galaxies, and these 
will also change the weak lensing spectrum (see [1^ for 
an application to one modified dark matter model). Fi- 
nally, the different evolution history of Aedm (a) may also 
have implications for the formation of nonlinear struc- 
ture. These further possibilities are beyond the scope of 
this work and will be further pursued elsewhere. 
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FIG. 4: (color online) The CMB spectrum of the model Ap = 
with the AEDM(a) parametrization given in Eq. (|43p and 
parameters a,; = 0.01 (blue curve) and ai = 0.05 (red curve) 
respectively. The black curve is the SCDM model. The other 
parameters are b = 1.001 and a/ = 1.0. 

B. The Case of Ap = 

In the case of Ap = 0, Eq. ((35l) simply becomes an 
algebraic equation for A^r: 

2 1 

-k Att = A^j^^ + HA^j^^ — -upEDMO, Aedm 

— —KphO^Ab — KprO^Ar- (44) 

For the parametrization of Aedm (a) we will still use 
Eq. In this case, because Eq. is algebraic, the 

only variable to propagate in time is wedm [cf. Eq. (j32p ]. 
and as in i^ lIV Al we could take wedm = prior to a.; and 
then evolve it using Eq. (|32p for a > a^. 

In Fig. 0] we have plotted the CMB spectrum for such 
a model with b = 1.001 and different choices of a^, from 
which we can see that the \ow-£ CMB power is very sen- 
sitive to both b and ai . The reason is that, although the 
right-hand side of Eq. (|44p is scale-independent thanks to 
our parametrization Eq. (|43)) . its left-hand side is scale- 
dependent through the fc^ factor. Consequently, on large 
scales (small k) the EDM anisotropic stress A^r could be 
very large and so the gravitational potential cf) is signifi- 
cantly different than that in standard CDM[c.f. Eq. ((25l) ]. 
The late ISW effect is then modified, which enhances the 
low-^ CMB power. On smaller scales, A^ is suppressed 
by fc^^ and its effects soon becomes negligible, explain- 
ing why the high-£ CMB power is not influenced. Also, 
the smaller Ui is, so the earlier the CMB evolution of the 
Ap — EDM model deviates from the standard CDM re- 
sult. Note that, in the case of standard CDM, the right 
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FIG. 5: The CMB power spectrum for the model in which 
Aedm (a) is parameterized as Eq. (|43p . both Ap and A^ 
are nonzero and satisfy a relation as described in the text. 
The solid, dashed and dotted curves represent the cases of 
b = 1.0, 1.1, 0.9 respectively, and they cannot be distinguished 
in the figure because they experience identical gravitational 
potential evolution. The other parameters are ai — 0.005 and 
af = 1. 

hand side of Eq. (|^^ vanishes identically, and so there is 
no influence on the last ISW. 

It is then clear that the Ap — EDM model with 
a general scale- independent parametrization of Aedm (a) 
is problematic and already stringently constrained. A 
possible way out of this trouble is to drop the scale- 
independence of Aedm (a) ■ We have numerically checked 
that if on large scales Aedm grows as Aqdm, then the 
low-i boosts of the CMB power as shown in Fig. 0] dis- 
appear and one recaptures the standard CDM results of 
ISW effect. Another possibility is to have both A^r and 
Aedm (a) scale- independent and also include Ap; this will 
be considered in § IIV CI 



C. General Ap and Att 

If the EDM arises from modifications to standard gen- 
eral relativity, then in general neither Ap nor A^r would 
be exactly zero. In this case there are many more pos- 
sibilities because Ap and A^ cannot be uniquely solved 
from Eq. ((35)) . Of course, the evolution of (j) might also 
be modified normally and this could be used to constrain 
Ap and A^. Indeed, we could use the freedom to choose 
Ap and A^ so that the evolution of cj) [cf. Eq. is 
exactly (or nearly) the same as that in standard CDM. 
Let us consider such an example now. 

In order that the gravitational potential evolves in the 
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FIG. 6: (color online) The quantities = Ap/Aedm (solid 
curves) and = A.^/Aedm (dashed curves) as functions of a. 
The black line is for SCDM (c^ = = 0), and the red/blue 
lines are for the model described in the text with h — 0.9, 1.1 
respectively. All curves are for the scale k — 0.01 Mpc~^. 



same way as in standard CDM, we require the contribu- 
tion on the right-hand side of Eq. (|36p from the EDM to 
vanish, which leads to the following evolution equation 
for 



4pr 



A" 

- m' - ^^^^ 



3pn 



4/9,. + 3p„ 



-HA' 



(45) 



ipr + 3p, 



Apr + 3pr, 



where pm and pr are respectively the energy densities 
of non-relativistic (baryons plus EDM) and relativistic 
(photons plus neutrinos) matter. Here note that the time 
evolution of A^^ is driven by Ap, which itself evolves in 
accord with Eq. psp . which is driven by A.^. plus the 
parametrized AEDM(a) (with its time derivatives). So, in 
the numerical calculation we have four more variables to 
propagate: Ap, dedm, A^r and A^. 

Note that if Ap vanishes identically and A^ — ini- 
tially then A.^. will remain zero all the time; also if A^ 
vanishes identically then so does Ap. These restrictions 
indicate that the two cases we have considered in above 
subsections cannot give rise to the same evolution of (f) 
and different evolution of Aedm to that predicted by 
standard CDM at the same time. The general case with 
Ap, A^ ^ 0, however, is able to do this, as is shown in 
Fig-El There, we plot the CMB power spectra of this gen- 
eral case with different values of b = 1.1, 1.0, 0.9 respec- 
tively, and they are totally indistinguishable because the 
evolutions of (f) are identical, leading to identical (zero) 
ISW effects. In this figure we have chosen = 0.005 and 



FIG. 7: (color online) The CMB spectrum for the general case 
[c.f. Eq. (|45|l ] with Aedm (a) parameterized as in Eq. (|43|l . 
Here we have chosen at — 0.0002 and af — 0.002. The blue 
dashed, black solid and red dotted curves represent the cases 
for b — 1.2, 1.0 and 0.8 respectively, as also shown beside the 
curves. The large-angle {£ < 100) CMB powers are indis- 
tinguishable for the curves as in Fig. [5] due to the identical 
evolutions of (j> and identical ISW effects. 



1 ' 1 ' 

CMB EE Spectrum 




FIG. 8: (color online) The CMB polarization spectrum for 
the model with the same parameters as in Fig. [T] The red 
dotted, black solid and blue dashed curves are the cases b = 
0.8, 1.0, 1.2 respectively. 
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a/ = 1 so that we only require the EDM to deviate from 
CDM after last scattering. It must be emphasized that, 
although the three curves are indistinguishable from each 
other, they do correspond to different EDM properties. 
To show this point clearly, we display in Fig. [6] the quan- 
tities = Ap/Aedm and = Att/Aedm, which char- 
acterize the importance of Ap and of A^ respectively, 
for the models. It is obvious that in the 5 7^ 1 models 
these quantities could be very different from the standard 
CDM value (0), and this can be understood intuitively as 
follows: the pressure perturbation term Ap acts against 
gravitational collapse, when 6 > 1, meaning that Aedm 
grows faster than Acdm early on and more slowly later; 
Ap needs to be negative early on and positive later (and 
vice versa). The anisotropic stress term A^r dissipates 
fluctuations in Aedm [c-f. Eq. (|35p ] and behaves simi- 
larly in the figure. 

We have thus seen that the observations of the CMB 
and matter power spectrum cannot rule out this general 
case (as long as is after last scattering) . Cross correlat- 
ing ISW with the galaxy distribution might help in this 
regard, but it is still limited for two reasons: first, it is 
contaminated by the effects of dark energy, and second, 
the cross-correlation data only exists for late times and 
cannot effectively constrain models like ours where the 
deviation from CDM occurs much earlier. 

If we choose the < 10~^, which means that the de- 
viation of EDM from CDM starts before last scattering, 
then although Eq. (|45p guarantees that the cj) evolution 
is not changed for arbitrary b, the CMB power spectrum 
will generally be different from that of standard CDM 
because the quantities <& and 5*, which are directly rele- 
vant for the primary CMB anisotropy, are modified since 



-0- 



In Fig. [7] we present such 



an example, for which we have used the parametrization 
Eq. (US]) of Aedm (a) with = 0.0002 and a/ = 0.002. 
This shows that in general the small-angle (high-^) CMB 
power will be different even though the large-angle power 
is fixed to be the same as in standard CDM. If the evo- 
lution of Aedm prior to last scattering is significantly 
different from that of CDM, then the deviation of CMB 
spectrum might be very large, and this is why the CMB 
could efficiently constrain such alternative gravitational 
dark matter theories as TeVeS and the f{K) models. 

Another possible constraint comes from the CMB po- 
larization. The linear polarization of CMB photons car- 
ries information about the gravitational potential be- 
cause it can only be generated from the quadrupole mo- 
ment of the CMB temperature perturbation, which is 
related to the anisotropic photon stress [tO, [zH , through 
Thomson scattering. Furthermore, at earlier times the 
rapid scattering of photons by electrons ruins the photon 
anisotropic stress (c.f. ^ IIII[) . consequently the polariza- 
tion only appear close to last scattering when significant 
photon quadrupole anisotropy can be produced. The lo- 
calization of the generation in time, and the dependence 
on the quadrupole moment, only mean that the polariza- 
tion directly refiects the gravitational potentials at the 



time of last scattering in a special way and so provides 
invaluable information about the EDM. This is in con- 
trast to the CMB temperature spectrum which depends 
on the gravitational potentials through the whole cosmic 
history up to now, and also on the photon monopole and 
dipole moments, which complicate the extraction of in- 
formation. In Fig. [5] we plot the CMB EE polarization 
spectrum for the same model as in Fig. [T] It can be seen 
that different parameters give quite distinct peak fea- 
tures. The CMB polarization was first detected in 2002 
72 1 , with the precision gradually improved since then 
73| . Although the precision at present is still insufficient 
to place stringent constraints on the EDM model, we ex- 
pect that future observations will change this situation. 

We could also slightly change Eq. ([45]) to make the 
evolution of (f> or ^E* the same as that for standard CDM, 
however clearly this is not possible for both $ and ^' be- 
cause in this case Att ^ 0. Because the primary CMB 
anisotropy is determined by both of these two variables, it 
will almost definitely differ between the EDM and CDM 
models. Furthermore, in the most general cases the evo- 
lution of (j) will be changed as well, which gives rise to 
different low-^ CMB spectra. Therefore, we expect that 
the CMB data will place stringent constraints on the gen- 
eral EDM models, and can be used to distinguish alter- 
native gravitational theories of dark matter. Conversely, 
any attempt of modification of general relativity which 
claims to have reproduced the observed large scale struc- 
ture must be confronted with the CMB fluctuation spec- 
trum and polarization to test its viability. 



V. SUMMARY AND DISCUSSION 

Motivated by the recent developments in producing the 
large-scale structure formation with alternative theories 
of gravity, we have considered the prospect of using the 
CMB to constrain such theories in new ways. We take 
the confrontation with the matter power spectrum data 
as a first test of the perturbed cosmological model in 
our alternative gravity theory and assume that this test 
has been passed. This is because the currently success- 
ful theories like TeVeS and f{K) really are only able to 
reproduce the observed LSS rather than all the pertur- 
bation observables, and such a simplifying assumption 
efficiently restricts our model space so that the degen- 
eracy problem is somewhat alleviated. As discussed in 
§ nil reproducing the observed LSS is a rather weak re- 
quirement on the theory because the LSS only reflects 
the density perturbations at late time. This means that 
there is still considerable freedom to choose the evolution 
history of the dark-matter density perturbation Aedm, 
which will generally lead to different predictions about 
other cosmological observables such as the CMB spectra 
and polarization. 

We discussed how this freedom in the evolution his- 
tory can be utilized by considering an example of the 
parametrization [c.f. Eq. (j43p ] of Aedm(i) which only 
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reduces to the familiar AcDM(a) result with the param- 
eter choice b = 1.0. As was shown in § [III the evolu- 
tion of Aedm is completely governed by the EDM stress 
history, which we quantify using the two variables Ap 
and Att. These same two variables also control the evo- 
lution of gravitational potential, which is important for 
the CMB power spectrum. This implies that the CMB 
is an ideal test bed for EDM models. In reality, we have 
one more freedom because once the Aedm (a) is specified. 
Eq. ([55]) becomes a single equation for the two variables 
Ap and A^r, so, in § IIV[ we considered three separate 
cases: (I) A^ =0, Ap ^ 0, (II) A^ ^ 0, Ap = and 
(III) Ajr 7^ 0, Ap 7^ 0. For simplicity, we have also made 
the following assumptions: (1) the evolution of Aedm 
is scale independent, as is implied by the observed LSS, 
(2) the cosmological constant (or other form of explicit 
dark energy violating the strong energy condition) is ne- 
glected to simplify the numerics, (3) the EDM equation of 
state parameter, weum, is assumed to be small enough so 
that we can assume the standard CDM background evo- 
lution holds; and, (4) only adiabatic initial conditions, 
a scale-independent 71^ = 1 primordial power spectrum, 
and effectively massless neutrinos are considered (we will 
return to these assumptions below). 

For case I, the standard CDM relation 'J = — $ = 
still holds but the evolution of these potentials could 
be changed. If the deviation of EDM from CDM only 
occurs after the last scattering, then this change in cj) 
mainly modifies the ISW effect and, hence, the low-^ 
CMB power. If the EDM starts to deviate from CDM 
before last scattering, however, then the early changes in 

and $ would modify the CMB acoustic peak features 
in a complicated manner, as we described qualitatively in 
§ mil and § IIV Al The evolution of Ap can be understood 
qualitatively: since the effect of the pressure perturba- 
tion is to counteract gravitational collapse, in order that 
Aedm grows faster than Acdm we need a negative Ap, 
and vice versa (c.f. Fig.E]). 

The case II evolution is problematic as long as we stick 
to the above simplifying assumption (1), i.e., requiring a 
scale- independent growth of Aedm, because in this case 
Eq. ([55l) implies a scale dependent A^ (fc^A^ is indepen- 
dent of k), which diverges on large scales (small k). As 
we have seen in § lIVBi this leads to a strong late-ISW 
effect, blowing up the low-^ CMB spectrum even if the 
EDM only differs slightly from standard CDM . Further- 
more, a natural EDM model with Ap = and A^r ^ 
does not exist in the literature, and this model is differ- 
ent from of case I, which arises naturally in f{K) models 
with ci3 = 0. 

Of course, the most general situation is our case III, 
where both Ap and A^ are nonzero. Needless to say, 
this generality also makes its exploration more difficult. 
One could, however, use the extra degree of freedom here 
in model constructions. If we fix both the evolution of 
Aedm (a) and that of the potential (p, then A^ and Ap 
can be solved uniquely at the same time. Such a con- 
struction has the advantage that we can fix the evolution 



of (j> to be identical to that in standard CDM so that the 
ISW effect is not changed at all. If the deviation of EDM 
starts after last scattering, then such a model is degener- 
ate with standard CDM even after CMB and LSS data 
are taken into account, and furthermore the CMB-galaxy 
cross correlation data is inefficient in constraining it. Al- 
lowing EDM to deviate before last scattering, however, 
will almost surely predict different CMB peak features, 
because in this case \1/ ^ — <(> ^ and the acoustic oscil- 
lations of photon-baryon fluid are changed with respect 
to standard CDM (c.f. i) IIII[) . Dropping the requirement 
on the evolution of cj) makes the situation even worse, be- 
cause in this case the ISW effect and low--^ structure of 
the CMB spectrum are also changed. We stress that the 
CMB polarization also provides invaluable information 
for constraining EDM models, because it depends only 
on the quadrupole moment of the photon temperature 
perturbation O (the photon anisotropic stress), and thus 
only on the physics close to the time of last scattering. 

Throughout this paper we have chosen to parametrize 
Aedm (a) as in Eq. (|^ . This is fairly simple and suf- 
ficient for our purposes as we only aim to show the im- 
portance of CMB data in constraining alternative grav- 
itational dark matter theories, but not to constrain any 
specified model or given parametrization exactly. There 
are good reasons why more detailed parametrizations 
are needed for more precise future studies. First, al- 
though Acdm oc a is a good approximation in the matter- 
dominated epoch, it breaks down when the contribution 
from radiation is still significant (including the time prior 
to and around last scattering) and when a cosmologi- 
cal constant is included. Second, because the physics 
at the time of matter-radiation equality is relevant for 
EDM models which deviate from CDM earlier, we some- 
times need to choose to be smaller than the matter- 
radiation equality Ooq, which means that the evolution 
of AEDM(a) or Acdm (a) cannot be completely scale in- 
dependent (remember that agq is relevant for the bend- 
ing of the matter power spectrum [t^). Thus, future 
precise calculations, particularly those relevant for weak 
lensing and the galaxy-ISW cross correlation, should use 
parametrizations which reduce exactly to Acdm (a) in the 
ACDM model in appropriate limits. 

Notice that our logic here is different from other related 
considerations of generalized dark matter (CDM) [sl, 
[tsI . In Ref. [t^ , for example, the author parametrizes the 
stress sector of the GDM with three quantities, Wg, c^^^ 
and c^j^, which respectively characterize the EOS, pres- 
sure perturbation, and anisotropic stress of the GDM. 
The GDM effects on the CMB and matter power spec- 
tra are then studied by assuming some specific values of 
these quantities. Our approach tackles the problem from 
a different direction and is therefore complementary to 
those earlier works. 

One interesting issue is the inclusion of hot dark mat- 
ter, which is needed by TeVeS itself. This can certainly 
be achieved by setting special values for Wg, c\jf and c^^^ 
as in [tHI . In our approach we simply treat the EDM as 
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a " black box" which is completely controlled by separate 
conservation equations. By parametrizing AEDM(a) we 
are able to extract information on the pressure pertur- 
bation Ap and anisotropic stress A^r from these conser- 
vation equations without knowing exactly what is in the 
box - it can be mixtures of modified gravitational effects 
and neutrinos, or something else entirely. 

The other issue is related to the initial conditions. In 
this work we mainly focus on the adiabatic initial con- 
ditions as in [7^. However, it is well known that the 
adiabatic initial condition is far from the only possibil- 
ity, and there can be four regular isocurvature modes [76| . 
These isocurvature modes are predicted by many theo- 
retical models and may even correlate with the adiabatic 
one. In contrast to the adiabatic mode, the isocurva- 
ture mode excites sinusoidal oscillations of the photon- 
baryon fluid and so predicts the first CMB acoustic peak 
to be around £ ~ 330 rather than 220; consequently 
a pure or dominating isocurvature initial condition is 
incompatible with basic CMB observations. Nonethe- 
less, a subdominant contribution from isocurvature ini- 
tial conditions actually degenerates with other cosmolog- 
ical parameters and cannot be ruled out by the current 
data [73, H H US- In our EDM model, if the EDM 
evolves differently from CDM early in the radiation era, 
then Aedm does not evolve adiabatically and the en- 
tropy perturbation 5edm — Aedm — j^-y is nonzero. In 
this case, one might naturally expect that there should 
be isocurvature modes in the initial conditions. A de- 
tailed calculation of the consequences of introducing such 
modes, however, generally requires thorough searches of 
the new parameter space (which is enlarged compared to 
the standard case because now the amplitude, tilt of the 
isocurvature modes and their correlation with the adi- 
abatic mode must be taken as free parameters) like in 
[73, H [zi, [13] and is beyond the scope of this work. 
Moreover, allowing different initial conditions (especially 
tilts of the isocurvature mode which are significantly dif- 
ferent from 1) can change the shape of matter power spec- 



trum, which means that the parametrization of Aedm (a) 
[c.f. Eq. should be scale dependent for the same rea- 
sons as discussed above. In this work we adopt a more 
conservative approach by assuming that the EDM evolves 
like standard CDM at earlier times and so adiabaticity 
is a natural choice. 

In conclusion, we propose ways to use the CMB in 
order to constrain those alternative gravity theories for 
dark matter which claim to be compatible with the ob- 
served LSS. We find that the CMB temperature and po- 
larization spectra are good discriminators between these 
theories in general, especially when they deviate from 
the CDM paradigm before last scattering. If the devia- 
tion starts after last scattering, however, there can exist 
EDM theories which are degenerate with respect to the 
standard ACDM model and cannot be distinguished by 
CMB and matter power spectra. This point is particu- 
larly interesting from the viewpoint of model construc- 
tions. Our results also indicate that the stress properties 
of dark matter, which determine the evolutions of both 
density perturbations and gravitational potential, can be 
studied and significantly constrained with existing and 
future data by using just the general conservation equa- 
tions and without specialising to any specific theoretical 
model. 
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